function y = fun(x)
    y=(-0.000007568286*x.^3 + 0.007775195238*x.^2 -2.450233333333*x+251.299047619048);
end